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Abstract 

We study the steady state motion of bubble walls in cosmological phase tran- 
sitions. Taking into account the boundary and continuity conditions for the fluid 
variables, we calculate numerically the wall velocity as a function of the nucleation 
temperature, the latent heat, and a friction parameter. We determine regions in the 
space of these parameters in which detonations and/or deflagrations are allowed. 
In order to apply the results to a physical case, we calculate these quantities in a 
specific model, which consists of an extension of the Standard Model with singlet 
scalar fields. We also obtain analytic approximations for the wall velocity, both in 
the case of deflagrations and of detonations. 

1 Introduction 

Phase transitions of the universe may give rise to a variety of cosmological relics, such as 
the baryon asymmetry of the universe p], cosmic magnetic fields |2J, topological defects 
[3], inhomogeneities [I], and gravitational waves [Hi El Ej. To be observable, several of these 
relics depend on the strength of the phase transition. In a first-order phase transition, 
bubbles of the stable phase nucleate and grow inside the supercooled phase. The expansion 
of bubbles provides a departure from thermal equilibrium, which is generally required to 
generate cosmological remnants. The mechanisms of relic generation depend in general 
on the motion of the bubble walls (either because they are based on charge transport near 
the bubble walls, on the collisions of the walls, or on the turbulence they produce). As 
a consequence, an important parameter in the generating mechanisms is the velocity of 
bubble expansion. For some relics, higher velocities give a stronger effect. This is the case, 
e.g., of gravitational waves. In other cases, on the contrary, the amplitude or abundance 
of the relic peaks at some value of the bubble wall velocity. This happens for instance in 
the case of electroweak baryogenesis. 
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A first-order phase transition requires the system to have two phases with different 
equation of state (EOS), which coexist in some temperature range. Thus, the high- 
temperature phase has a pressure p + (T), and the low-temperature phase has a pressure 
P-(T). At the critical temperature T c , the two phases have the same pressure but different 
entropy and energy. The entropy density is given by s(T) = p'{T), where a prime denotes 
derivation with respect to T, and the energy density is given by p(T) = Ts(T) —p(T). 
The latent heat L is defined as the difference between the energy densities of the two 
phases at T = T c . Therefore, we have p+(T c ) = P-(T C ) and L = T c (p' + (T c ) —p'_(T c )). 
Furthermore, all these quantities can be derived from the free energy density T by the 
relation p(T) = — JF(T). 

In a finite temperature field theory, we have in general a scalar field <f), the Higgs field, 
which plays the role of an order parameter for the phase transition. We shall assume 
this is the case, although most of our conclusions will be independent of this assumption. 
The free energy density is the finite temperature effective potential, which is a function 
of and T. In some range of temperatures, the free energy has two minima separated 
by a barrier. The value of <fi in each minimum defines the two phases. We shall assume 
that the high-temperature minimum is <fi = 0, and the low temperature minimum has a 
value 4>m{T) 7^ 0, which is the case for a symmetry breaking phase transition. A phase 
transition with a large value of the order parameter, <f> m (T) > T is usually called a strongly 
first-order phase transition. Thus, we have F+iT) = JF(0,T) and JF_(T) = jF(0 m (T),T). 
Below the critical temperature, the transition from (p = to <p — <Pm{T) occurs via the 
nucleation of bubbles [H A bubble is a configuration <p(x) of the field. We shall make 
the usual assumption that this configuration is a thin-walled sphere. A nucleated bubble 
grows because P-(T) > p+(T) for T < T c . 

The propagation of bubble walls in cosmological phase transitions has been extensively 
investigated Hd HH H21 H31 QH H51 HSl HZl HBl HSl I^S I21l I2H HSl I2H I2S]- In general, the 
stronger the phase transition, the larger the amount of supercooling. Therefore, one 
expects that the velocity of bubble walls will be larger for stronger phase transitions. 
Indeed, the pressure difference which pushes the wall of a bubble, Ap(T) = p~(T) —p + (T), 
grows as the temperature decreases. On the other hand, the velocity depends also on the 
friction of the wall with the medium and the latent heat that is released at the phase 
transition front. The friction is related to the microphysics, i.e., to the interactions of 
particles of the medium with the Higgs field in the configuration of the wall, whereas the 
latent heat is involved in the hydrodynamics, i.e., the bulk motions and reheating of the 
fluid near the wall. 

According to hydrodynamics, there are two kinds of steady state solutions for the 
propagation of the wall, namely, detonations and deflagrations. Detonations are super- 
sonic while deflagrations are generally subsonic. Each solution has different boundary 
conditions, and both may appear in a given model and may even coexist in some range 
of parameters. A numerical calculation of the evolution of the field 4>(x,t) [TO], suggests 
that the bubbles grow typically as subsonic deflagrations, although, in some regions of 
parameter space, the solutions switch from deflagrations to detonations, as the friction is 
decreased. In view of the cosmological consequences of the phase transition, it is impor- 
tant to determine the range of parameters where each solution may exist. Hydrodynamic 
considerations give general constraints on the allowed regions of thermodynamic param- 



2 



eters, but the existence of deflagrations and detonations depends strongly on the value of 
the friction. 

The literature on the subject of the bubble wall velocity can be roughly divided in two 
groups: those papers which calculate the friction for a given model, assuming nonrelativis- 
tic deflagrations and ignoring hydrodynamics [HI [121 [131 CHI CHE EE E] > an d those which 
consider hydrodynamics, but they either disregard microphysics or include the friction 
as a free parameter pUl HHJ, [Hi EH EU [22j [23j [24]. As a consequence, the analytical ap- 
proximations that are most commonly used in applications, were derived ignoring either 
microphysics or hydrodynamics. Indeed, in the context of baryogenesis the wall velocity 
is assumed to be of the form v w ~ Ap(T)/rj, where rj is a friction coefficient. This is a 
nonrelativistic approximation (thus appropriate for the case of deflagrations) which ne- 
glects hydrodynamics. In the context of gravitational waves, on the contrary, the wall is 
assumed to propagate as a detonation, and the so called Chapman- Jouguet hypothesis is 
further assumed, leading to a simple expression for the velocity, which depends only on 
the ratio between the latent heat and the thermal energy density of the plasma. 

The aim of this paper is to calculate the wall velocity in both cases, detonations and 
deflagrations, taking the friction into account. This will allow us to investigate under 
which conditions the various hydrodynamic solutions can exist, and to explore how these 
conditions are attained in a physical model. For these purposes, we shall first consider 
the equations for the discontinuity of the fluid variables across the wall, together with a 
simple parametrization of the friction force (which consists in introducing a damping term 
in the equation for the field <p) and a simple model (the bag equation of state) for the 
phase transition. These simplifications permit to write down a complete set of equations, 
which can be solved for the wall velocity as a function of three parameters, namely, 
a c = L/(4p + (T c )), where L is the latent heat and p+(T) is the thermal energy density of 
the high-temperature phase, ocn — L/(4p + (T/v)), where T/v < T c is the temperature at 
which bubbles nucleate, and rj/L, where rj is the friction coefficient. 

Solving numerically these equations, we shall determine the regions in the space of the 
aforementioned parameters for which deflagrations and detonations may or may not exist. 
In order to apply the results to a physical model, we choose as an example the electroweak 
phase transition in an extension of the Standard Model (SM) with scalar singlets, which 
provides a simple variation of the strength of the phase transition as a function of the 
fundamental parameters. For this specific model, we compute the quantities a c , a^, and 
rj/L and we calculate the wall velocity. To estimate the value of the parameter rj, we shall 
relate our parametrization of the friction to the microphysics calculations that exist in the 
literature. For the nucleation temperature T/v, we solve numerically the bounce action 
to obtain the nucleation rate. We shall also find analytical approximations for both 
propagation modes. In particular, we give a simple approximation for the detonation 
solution, which depends on the friction. 

The paper is organized as follows. In section [2] we write down the equations for the wall 
velocity, the fluid variables, and the friction for both propagation modes. In section [3] we 
show the numerical results for the wall velocity as a function of the relevant parameters, 
and we determine regions in parameter space where each kind of solution may or may not 
exist. In section H] we consider a specific model, for which we compute the values of the 
latent heat L, the friction rj, and the nucleation temperature T N , and we calculate the 
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detonation and deflagration velocities (when they exist) as functions of the fundamental 
parameters of the model. Finally, in section [5] we consider analytical approximations for 
the wall velocity. Our conclusions are summarized in section EJ 



2 The wall velocity 

To take into account hydrodynamics and microphysics in the calculation of the wall ve- 
locity, three basic ingredients are needed, namely, the continuity conditions that relate 
the fluid variables across the wall discontinuity, an equation of state, which relates the 
various fluid variables on each side of the phase transition front, and an equation for the 
effect of friction on the motion of the wall. 



2.1 Hydrodynamics 

In the rest frame of the phase-transition front, we have an incoming flow in the high- 
temperature phase, with velocity v + , and an outgoing flow in the low-temperature phase, 
with velocity t>_ (see Fig. [T]). The motion of the wall causes the temperatures T + and T_ 
on each side to be different. The continuity conditions for energy and momentum fluxes 
give two equations |26j, 



W- 



2 2 i 

-7_w_ + P-, 



(1) 



where 7 = — v 2 and w = p + p = Ts is the enthalpy density. Equivalently, we have 



(P+-P-) v + (p-+p+) 



P- 



(P++P-) 



(2) 



which can be readily solved for the velocities in terms of the thermodynamical quantities. 
Notice that these two equations have four unknowns, namely, the velocities v± and the 
temperatures T± on both sides of the wall. All other thermodynamical quantities are 
determined by the equation of state. In principle, the temperature T + of the supercooled 
phase can be calculated using the nucleation rate. To determine completely the system, 
one more equation is needed. This can be obtained from a microscopic calculation of the 
friction of the wall with the plasma. 



detonation 



deflagration 



Figure 1: The velocities on each side of the wall, in the reference frame of the wall. 

The above equations admit two kinds of solutions, one of them with < |t> + |, 
called detonation, and the other with |i>_| > |t> + |, called deflagration (Fig. [p. This can 
be seen by using a simple model in which the low-temperature EOS is that of radiation, 
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P- = P-/3, and the high-temperature EOS is that of radiation plus vacuum energy, i.e., 
p+ = p + + e, with p + = p+/3 — e. From Eqs. (JIJ) or (j2J) one obtains [T8] 



-±- + — ± 



— + — 

6«_ ~ 2 



a 2 + ia 



1 



a 



(3) 



where a = e/p + . The solutions with the plus and the minus signs correspond to det- 
onations and to deflagrations, respectively. We plot both in Fig. [2] for a given value 
of a. It turns out that, for a detonation, we have \v + \ > c s+ , and for a deflagration, 
\v+\ < c s+ , where c s ± is the speed of sound on each side of the wall. For the above simple 
equations of state, the speed of sound in both phases is that of a relativistic plasma, 
c s ± = (dp±/dp±) 1/2 = l/y/3. 
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Figure 2: v + vs t>_ according to Eq. (j3J), for a = 0.1 



In the reference frame where the matter far in front of the bubble wall is at rest and 
the wall moves, the matter at the center of the bubble (far behind the wall) must be also 
at rest. As a consequence, a single front (the phase transition front) is not sufficient to 
satisfy all the boundary conditions. As a result, it turns out that a deflagration front must 
be preceded by a shock front (Fig. [31 left). Between the shock and the deflagration fronts, 
the matter is compressed and has a finite velocity. Outside this region the matter is at 
rest. On the other hand, a detonation front hits fluid which is at rest, the density increases 
suddenly, and the phase transition front is followed by a rarefaction wave (Fig. [31 right). 
Comparing with the reference frame of the wall, we see that the deflagration wall velocity 
is given by v w = |t>_|, while the detonation wall velocity is given by v w = \v+\ . According 
to the inequalities above, a detonation is always supersonic, v w > c s+ . A deflagration may 
in principle be supersonic or subsonic (it is subsonic with respect to the fluid in front of 
it, which is moving in the direction of the wall). 
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Figure 3: Schematically, the profile of the density or the velocity of the fluid, induced by 
the moving wall. 



2.1.1 Strong and weak solutions 

Depending on the velocity V- of the outflow, deflagrations are divided into weak (|f_| < 
c s _), Jouguet (\v-\ = c s _) and strong (|t>_| > c s _) deflagrations. Thus, strong deflagra- 
tions are supersonic, and weak deflagrations are subsonic. Detonations are also classified 
into strong (|d_| < c s _), Jouguet (|t>_| = c s _), and weak (|t>_| > c s _) detonations. Since 
the incoming flow is supersonic, strong detonations perturb the fluid more strongly than 
weak detonations. 

The various kinds of solutions have been extensively investigated [101 [HJ [191 EH EE] . 
Using thermodynamical arguments, one can find regions in parameter space where each 
solution is allowed. In particular, the non-negativity of entropy production may rule out 
the whole family of detonations or deflagrations [21 J . However, such arguments do not 
rule out a particular type (i.e., weak, Jouguet, or strong) of detonation or deflagration 
|22j . In the case of deflagrations, in Ref. [22J it is argued that strong deflagrations are 
mechanically unstable, although in principle there is no physical reason against their 
existence. A special type of supersonic Jouguet deflagration, which is followed by a 
rarefaction wave, was shown to exist in Ref. [23J. 

In the case of detonations, it is known [2S] that, for the case of chemical burning, 
detonation bubbles can only grow as Jouguet detonations. Remarkably, this condition 
(which is referred to as the Chapman- Jouguet hypothesis) means that a microscopic 
calculation of the wall velocity is irrelevant, since the system is completely determined 
by the energy-momentum conservation and the boundary conditions. In Ref. [IB] it was 
argued that the Chapman- Jouguet hypothesis should be valid also for cosmological phase 
transitions. Setting t>_ = c s in Eq. ([3]) leads to a very simple formula v w = vj(a) for 
the wall velocity. However, it was shown in Ref. [22] that, in contrast to the case of 
chemical burning, the Chapman- Jouguet hypothesis does not hold in the case of phase 
transitions. It was also shown in Refs. [T5[ [22] that strong detonations are impossible, 
since they cannot satisfy the boundary condition of a vanishing fluid velocity at some point 
behind the bubble wall. Nevertheless, weak detonations are possible, and the velocity 
will be determined in general by microscopical processes. Thus, in a realistic model the 
detonation velocity will depend on a friction parameter rj as well as on the hydrodynamic 
parameter a. The wall velocity will always satisfy the condition v w (a, 77) > vj (a) (as can 
be seen in Fig. [2]). Hence, v w may give the Jouguet result only in some limiting cases. 
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2.1.2 Detonations vs deflagrations 



It is well known that, if the supercooling is not considerable and the friction is large 
enough, the bubbles will grow as deflagrations, whereas detonations demand small values 
of the friction and more extreme conditions (i.e., strong supercooling) [TOl [T9| |2T]. It is 
often argued, though, that the propagation of the deflagration wall becomes turbulent 
and eventually turns into a detonation due to shape instabilities of the wall. As discussed 
in Ref. [27] for the case of the QCD phase transition and in Ref. [28] for the case of the 
electroweak phase transition, small wavelength perturbations of the deflagration front are 
stabilized by the interface tension, while large ones grow exponentially. However, a more 
detailed study [29] shows that the strong dependence of the wall velocity on temperature 
actually stabilizes the large scale perturbations. The stability of the deflagration front 
thus depends on several quantities (these are essentially the latent heat, the value of the 
wall velocity and its dependence on temperature). It was found, in particular, that the 
deflagration front is stable under perturbations at all scales for velocities above a certain 
critical value v c . For instance, in the (unrealistic) case of the minimal Standard Model 
with a Higgs mass m# = 40GeV, the critical velocity is v c ~ 0.07 [29]. Below this velocity, 
perturbations are unstable on scales larger than a critical scale A c . However, depending on 
the parameters, the value of A c may be much larger than the size of bubbles. Furthermore, 
the characteristic time for the growth of the instabilities may be larger than the duration 
of the phase transition. Therefore, in the general case, the deflagration front will be 
probably stable. A study of the hydrodynamical stability of bubble growth is beyond the 
scope of this paper. 



2.2 Microphysics 

In any case, it is clear that the determination of the wall velocity requires additional input. 
To consider the effect of friction together with the fluid equations ([TJ) , it is customary to 
introduce a damping term of the typical form u^d^ip in the equation of motion for the 
Higgs field rjO], 

<W + 9:F ^ T) + (T c fj) u'drf = 0. (4) 

Here, the dimensionless damping coefficient fj is a free parameter. Nevertheless, we can 
relate fj to the actual friction coefficient obtained from microscopical calculations in a 
specific model, as we shall see in section 14.4.11 For simplicity, we shall consider a planar 
wall moving along the z axis. In the frame of the wall we have 

- 0" + dF/dcj) + T c f) U( j)' = 0, (5) 

where u = jv, with v the fluid velocity in the z direction. 

Multiplying by <$' (z) and integrating over — oo < z < +oo, the first term disappears. 
Then, if we neglect the variation of temperature with z, the fluid velocity is v — —v w , and 
we obtain 

T (0, T)-T (4> m , T) = T c fjajv w , (6) 

where 

a = I V (z)} 2 dz (7) 
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is the surface tension of the bubble wall. Equation ([6]) gives the well known nonrelativistic 
approximation for the wall velocity 



v w nAF(T)/r), (8) 

where 

77 = rjT c a, (9) 

and AjF is the free energy density difference JF + (T) — JF_ (T), which gives the pressure 
difference Ap = p^ (T) — p + (T) between the two phases. 

If we now take into account hydrodynamics, we can still integrate Eq. (jSJ), but we can- 
not ignore the fact that the temperature and velocity of the fluid depend on z. Before mul- 
tiplying by 4>' and integrating, we use the equation dJ^/dcj) = dF/dcfi — (dT/dT) (dT/d(f>) 
[TO] . Then, proceeding as before, we obtain 

T (0, T+) - T (0 m (T_) , T_) - / ^dT + f,T c \ \<p' {z)f u (z) dz = 0. (10) 



T- 



The first two terms are just JF + (T + ) — T- (T_) = p- — p + . To simplify the calculations, 
some approximations are necessary. In the first integral, the function dT j dT (0 (x) ,T (x)) 
varies from 

^ (0 m (T_) , T_) = ^ (T_) = - a _ (T_) (11) 

to 

^(0,T + ) = ^(T + ) = - S+ (T + ) (12) 

(we have used the fact that m and are minima of J 7 ). If (x) and T (x) change smoothly, 
we expect that a linear function will give a good approximation to the integrand dT /dT 
inside the thin wall. Therefore, the integral yields — (s_ + s + ) (T + — T_) /2. For the 
second integral, it is a good approximation to assume that 0' 2 is symmetric around the 
center of the wall at z = 0; hence, this can be approximated by a (u+ + u_) /2 pjj]. Hence, 
we have 

P+ - P- ~ \ (s+ + O (T+ - r_) + I (i; +7+ + U-7-) = 0, (13) 

where rj = fjT c a as before. 

Notice that the approximations we have made to obtain Eq. ( TTSl do not involve any 
assumption on the wall velocity nor on the equation of state. Hence, this equation can 
be used for the treatment of relativistic as well as nonrelativistic velocities, and is model- 
independent. Provided the various thermodynamical quantities are related by means of 
an equation of state, Eqs. (JTJ) and ( TlBl can be solved to obtain the velocities t>_ and 
v + as functions of the temperature T + in front of the wall. The result will depend on 
parameters of the theory that appear in the EOS, such as the latent heat, and on the 
friction coefficient. In practice, however, it is very difficult to solve these equations unless 
we use a simple model. 
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2.3 The equation of state 

As we have seen, the relation (jHJ) between v + and t>_ can be derived from Eqs. ([I]), using 
the relations 

V- = P-/3 (14) 

for the low-T phase, and 

P+ = P+ + £, P+ = P+/3 -e (15) 
for the high-T phase. In order to obtain a second relation between v + and v _ from Eq. 
( fT3|) . the EOS for each phase has to be more specific. The thermal energy densities p_ (T) 
and p + (T) must be of the form aT 4 , where a = 7r 2 g/30 and g is the number of relativistic 
degrees of freedom. Hence, the simplest possibility is to assume that both phases have 
the same a. Such a simple model was considered in Ref. [18J. This model, however, fails 
to describe a realistic phase transition if the vacuum energy e is a constant. Indeed, the 
critical temperature is defined as that T at which the pressures of the two phases are 
equal, p + (T c ) — p_ (T c ) = 0. However, with p + = aT 4 /3 — e and p_ = aT 4 /3, we have 
p + (T) — p_ (T) = — e at any T. This would be possible if e were a function of T. With e 
constant, this model does not have a critical temperature. 

The existence of a critical temperature T c , below which bubbles nucleate and grow, 
constitutes an important distinction between the case of a phase transition and that of 
chemical burning. In a phase transition, supercooling is needed for bubbles to grow. 
In the first place, the temperature must decrease sufficiently below T c so that bubbles 
nucleate. In the second place, even if for some reason (e.g., by means of inhomogeneous 
nucleation) bubbles nucleate at T = T c , the condition T < T c is still necessary to achieve 
a pressure difference p~(T) > p + (T), so that bubbles can grow. In the case of chemical 
burning, the speed of the reaction increases with temperature, and the temperature must 
rise sufficiently for the combustion to proceed. If the reaction is strongly exothermic, it 
is sufficient to rise the temperature at a single point. The heat that is released by the 
reaction at that point rises the temperature of the surrounding gas, and the reaction may 
extend over the whole gas [26J. As a consequence, the higher the released energy, the 
larger the velocity of the combustion front. In the case of a phase transition, the release 
of latent heat is not necessary for the process to continue. On the contrary, the reheating 
of the supercooled gas slows down the phase transition, as the temperature approaches 

v;0 

To keep the EOS as simple as possible, we shall consider e constant, and different 
values a + and a_ in the two phases. This gives the well known bag equation of state, 
which is appropriate to a system with a first-order phase transition: 

p + = a+T 4 + £, p + = a + T 4 /3 — £, (16) 
p_ = a_T 4 , p_ = aJT /3, 

and the entropy is given by s± = |a±T 3 . Since Eq. ([3]) was derived using the relations 
( I14|ll5p . the result is still valid as a function of a = s/p + , with p + = a + T 4 . However, 

lr The process we have just described corresponds to a deflagration or slow combustion. For the 
detonation case one can argue similarly. Essentially, in the case of chemical burning, the combustion 
takes place behind the detonation front because the temperature is risen by the front. Evidently, a phase 
transition does not need such a temperature rise. 
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the parameter e has a different meaning in each model. Indeed, with a + = a_ the energy 
released in the phase transition is just e. In a realistic model the vacuum energy e does 
not coincide with the latent heat L = p + (T c ) — p_ (T c ), since entropy is released together 
with vacuum energy. For the bag equation of state, the critical temperature is readily 
obtained by equating the pressures p + and p_. We have 

(a+ -a^)T c 4 = 3e, (17) 

from which we obtain the latent heat 

L = As. (18) 

In fact, in a realistic model, the vacuum energy density does not vanish immediately after 
the phase transition. Hence, we will have in general two parameters e±, one for each 
phase. Nevertheless, this only amounts to replacing e = e + — e_ in Eqs. (117|I18I) and in 
the definition of a. The expression of a in terms of the latent heat, a = L/(4p+), remains 
unchanged. For simplicity, we shall set e_ = in the following. 

Notice that the positivity of the pressure at T = T c (i.e., p + = p_ > 0) implies that 
the ratio 

a c = e/(a + T?) (19) 

cannot be arbitrarily large, namely, a c < 1/3. This relation between the vacuum and the 
thermal energy densities is valid beyond this simple model [30] . 

Using the fluid conditions (DO) and the equations of state ( fT6l) . the friction equation 
( [TBI becomes 

with 

P+-P- ^ v + v - / -n 

— = = a+- , 21 

p + 1 — 6v + v - 

/rn \ 3 rri / \ 1/4 

s_ a_ / 1 \ i_ f a + p. 



, 7 , ■ r , i (22) 

and 

— = l-a +T7 r , 23 

where we have used instead of a the notation 



a 



+ ^e/{a + Tl). (24) 



The velocities are further related by Eq. ([3]), with the + sign for detonations and the - 
sign for deflagrations, and the ratio a_/a + is a parameter of the model. Using Eqs. ( fTTl) 
and f|T9|) we can write 

a_/a + = 1 — 3a c . (25) 

It can be seen from Eqs. (T2J) that v + v_ < 1/3. Hence, according to Eq. ( J2T1) . we have 
p + (T + ) — p_ (T_) > 0. This means that the temperature difference that is established 
around the wall inverts locally the relation p_ (T) > p + (T) , which holds for T < T c . 
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We can solve the above equations for v+ or t>_ as a function of T + (or, equivalently, 
aq_) and the parameters r\ and a c . In the case of detonations, the temperature T + in front 
of the bubble wall is just the temperature Tjv at which bubbles nucleate, so we have 

a + = a N = — (detonations), (26) 
a + l N 

and the fluid on that side of the wall is at rest, so 

v w = \v+\. (27) 



Thus, we only need to solve Eqs. (I20H23P to obtain v w . This makes this case in general 
simpler than the deflagration case^l 

In the case of deflagrations, we must relate the temperature T + of the fluid in the 
reheated region between the phase-transition front and the shock front, to the nucleation 
temperature T N of the fluid beyond the shock front (see Fig. [3]). To do that, we must 
consider the fluid conditions ([T]) for the shock front discontinuity. These are simpler than 
those for the wall, since we have the same phase on both sides of this front. Calling v\ 
the velocity of the fluid in front of the shock and t>2 the velocity behind it, we have, in 
the reference frame of the shock front |26|, 



1 ( 3Ti + Tfc 



1/2 



M = 7=3 UtjTt? ) • H = 3RT (28) 

In the "laboratory" frame, the fluid is at rest behind the phase transition discontinuity, 
and also ahead of the shock front, so we have 

v w = \v.\, (29) 

whereas v\ gives the velocity of the shock front, V\ = — f s h- The velocity of the fluid 
between the two fronts can be computed from the velocities v + and i>_, or from V\ and 
v 2 . Equating both results we find the relation between T + and T N , 

^ V+ ~ V - (deflagrations). (30) 



^/{3a N + a+) (3a+ + a N ) 1 - v+V- 

Thus, we can find v w (a c , a^, rj), both in the case of a detonation and of a deflagration. 
Notice that the parameter a c gives the ratio of latent heat to thermal energy density 
at T = T c . Thus, for a c fixed, the parameter gives a measure of the amount of 
supercooling, since a^ja c = {T c /Tn) A . We have seen that < a c < 1/3. On the other 
hand, is only bounded below, > a c . Nevertheless, in general we will have T N ~ T c 
and a at < 1. We could have ^> a c only in a phase transition with extremely strong 
supercooling. 



2 If we wanted to find the temperature or velocity profiles inside the bubble, then we would have to 
consider the equations for the fluid behind the wall and impose boundary conditions at the center of the 
bubble. 



11 



3 Numerical results 



We have solved numerically Eqs. (I20H23P and (jBJ), with Eqs. (I26|27p for detonations and 
Eqs. ( I29f3"0~l) in the case of deflagrations. In this section we present the numerical results, 
considering independent variations of the parameters a c , and rj. Nevertheless, we 
shall make variations around two reference sets of parameters, corresponding to different 
points in the physical model of section 0] (indicated by crosses in Fig. [TTi) . We choose 
a c pa 1.27 x 1CT 2 , «at pa 1.48 x 1CT 2 , rj/L pa 4.03 x 1CT 2 , which are obtained when the 
fundamental parameters of the model are such that the phase transition is rather strongly 
first-order, i.e., the order parameter is (p^/T N pa 2; and a c pa 3.19 x 10~ 3 , pa 3.25 x 1CT 3 , 
rj/L pa 2.48 x 10~ 2 , obtained from a weaker phase transition, but still with (p^/T^ pa 1. 



3.1 Detonations 

Let us consider a strong phase transition (4>n/Tn ~ 2). In Fig. @] (left panel), we have 
fixed the values of a c and rj/L, and we have varied a^, starting from = a c , i.e., we have 
decreased the supercooling temperature T N below T c . For close to a c , no detonation 
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Figure 4: Detonation solutions for a c pa 1.27 x 10 -2 . Left: the velocities as functions of 
a N for rj/L pa 4.03 x 10~ 2 . Right: the velocities as functions of rj/L for a N pa 1.48 x 10~ 2 . 
Solid lines indicate the velocity v w of the wall, dashed lines the velocity v _ of the outgoing 
fluid, and dotted lines the speed of sound (lower curves) and the Jouguet velocity (upper 



curves . 



solutions exist (i.e., some supercooling is required). When is large enough we find two 
solutions for the detonation wall velocity (solid lines). One of them (lower solid curve) is 
close to the Jouguet velocity Vj (ckjv), which is indicated by the upper dotted line. This 
solution, however, does not seem to be physical. In the first place, the velocity decreases 
with the supercooling until it reaches the Jouguet point. Beyond this point, the velocity 
increases, but the solution has become a strong detonation, since the velocity t>_ (lower 
dashed curve) crosses below the speed of sound (indicated by the lower dotted line). The 
other solution, in contrast, is always a weak detonation (the corresponding velocity V- is 
indicated by the upper dashed curve). As a matter of fact, this detonation is quite weak, 
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as the outflow velocity t>_ is very close to the inflow velocity t> + = v w , which implies that 
the fluid is almost unperturbed by the wall. 

This behavior can also be seen if we fix a c and cxn, and plot the wall velocity as a 
function of r\ (Fig. HJ right panel). If the friction is strong enough, no detonation solution 
exists. For lower values of the friction parameter, we have two solutions, but the lower 
branch corresponds to a velocity which decreases as the friction is decreased. Depending 
on the values of the parameters, this solution eventually becomes a strong detonation for 
low values of rj. In any case, as can be seen in the figure, the lower branch solution is 
much stronger than that of the upper branch. For these reasons, we shall discard this 
solution. 

3.2 Deflagrations 

Solving numerically Eqs. (l20H23p . (j2D, and (1501) . we find again two solutions. However, 
one of them must be discarded, since it has a wrong behavior as a function of the pa- 
rameters: the velocity decreases with the supercooling and increases with the friction, 
and is in general supersonic. The correct solution, on the other hand, goes to zero for 
small supercooling and strong friction. As an illustration, we plot both solutions in Fig. 
Elfor the case of a strong phase transition, as functions of a at and rj. In contrast to the 
detonation case, the deflagration solution always exists in the limit of small supercooling, 
«iv — ► oi c . On the other hand, the deflagration solution may not exist when the supercool- 
ing becomes too strong or when the friction becomes too low, as this example illustrates. 




Figure 5: Deflagration solutions for the same case as Fig. HI Solid lines indicate the 
velocity v w of the wall, dashed lines the velocity v + of the incoming fluid, and dotted lines 
the speed of sound. 

In the high supercooling and low friction range, we have the detonation solution. We 
have plotted the physical detonations and deflagrations together in Fig. El Solid lines 
correspond to the case of Figs. H] and El and dashed lines correspond to variation of 
parameters beginning from a weaker phase transition (with <Pn/Tn ~ 1). As can be seen 
in the right panel, for a c and fixed, a stronger phase transition gives, for the same 
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friction, a larger velocity (solid line). On the contrary, if we vary a/v (left panel), the solid 
line gives lower velocities. This is because in this plot we are considering the the same 
amount of supercooling for the two models. Therefore, the model with a larger value of 
ot c (i.e., of L) gives a lower velocity. Indeed, a larger L causes a larger reheating, thus 
slowing down the phase transition (both in the deflagration and in the detonation case). 




Figure 6: Deflagration and detonation solutions as functions of (left) and of r]/L 
(right) for two families of phase transitions. 

In Fig. [6] we see that different situations may arise as the parameters are varied. 
There are ranges in which detonations and deflagrations coexist, and there are values of 
«at and r\ for which neither solution exists. Fixing a c , we have calculated the maximum 
friction for which detonations exist, and the minimum friction that permits the existence 
of deflagrations. We have plotted these limiting values in Fig. [7] (solid and dashed lines, 
respectively). Since the curves cross, the plane is divided into four regions, in which there 
exist either deflagrations, detonations, both, or none. In the latter case, no stationary 
state of bubble expansion is reached. In the case in which both solutions exist, which of 
them is realized will depend on the initial conditions and on the stability of each solution. 
Figure [8] shows the four regions for different values of a c . 

It is difficult to determine in general which of these situations may arise in a physical 
case. For that aim we would need some relation between the parameters rj/L, a^, and 
a c . In the next section we explore the physical case by considering a particular model, for 
which we calculate the relevant parameters 77, L, and T/v as functions of the fundamental 
parameters of the model. Nevertheless, some general conclusions can be drawn from Fig. 
Notice that the key parameter for the existence of detonations is the supercooling 
parameter a^ja c : even for r\ = 0, detonations exist only if a^ja c is large enough. On 
the contrary, deflagrations exist up to a certain amount of supercooling, even for 77 = 0. 
Therefore, in a phase transition with small supercooling, we will have only deflagrations, 
no matter how small the friction may be. This is because hydrodynamics brakes the 
propagation of the wall, thus acting effectively as a friction. This can also be seen in 
Fig. El by the fact that, as we increase a c (i.e., L), the deflagration region grows and the 
detonation region reduces. 
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Figure 7: The regions of the parameters r\ and a at, for a c = 0.05, where detonations and 
deflagrations exist. Solid line: maximum friction for detonations. Dashed line: minimum 
friction for deflagrations. 

Observing Figs. H or [8] we see that, setting rj = we find, for each a c , the minimum 
otN for which detonations may exist and the minimum value for which deflagrations may 
not exist. Besides, the value at which the detonation and deflagration curves cross defines 
the minimum for which there could be no solution. Fig. [9] shows these three values 
of «iv as functions of a c . These curves divide the (a c ,a^[) plane in four regions. In the 
lower region, detonations do not exist, independently of the value of rj; in the next region 
detonations will exist if the value of i] is small enough; in the third region, deflagrations 
will not exist if 77 is too small; in the upper region, neither solution will exist in a certain 
range of the friction parameter. Deflagrations always exist in the two lower regions and, 
for large enough friction, in all the (a c , a at) plane. 

4 The wall velocity in a physical model 

So far we have calculated the wall velocity as a function of the latent heat L, the super- 
cooling temperature T N , and the friction coefficient rj. In a physical model these quantities 
are not independent. However, the simple EOS f|T6|) does not allow to calculate Tjy nor V- 
We shall now consider a model which permits to obtain these quantities as functions of 
a single fundamental parameter. The model we seek should be as simple as possible, and 
allow for phase transitions of different strength as the parameters are varied. It is well 
known that the electroweak phase transition is only a smooth crossover in the minimal 
Standard Model (SM), whereas many extensions of the model give a strongly first-order 
phase transition. Therefore, we choose the simplest extension of the SM, namely, adding 
to the model a gauge singlet scalar S with zero vacuum expectation value [31j [32j [33] . 
Besides the coupling to the Higgs, which is of the form 2h 2 S> S H , the field S may 
have a mass term fi 2 S^S and a quartic term \ s (S^S) . We will ignore the possibility 
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Figure 8: The four regions in the plane (ajv,^) (as in figure [7]) for a c = 10 3 (solid), 10 2 
(dashed), 5 x 10~ 2 (dashed-dotted), and 10 _1 (dotted). 



that cubic terms exist in the tree-level potential, which may produce a barrier between 
minima at T = 0. Recently, an extension of the SM with several real singlets Si has been 
considered [31] ■ These bosons constitute a hidden sector which couples only to the SM 
Higgs doublet through a term h 2 WHY^ Sf. If the fields Si do not have mass terms, so 
that they only get a mass from electroweak breaking, the phase transition can be made 
exceedingly strong. 



4.1 The effective potential 

Our model will consist of a scalar field <p (the background Higgs field, (H°) = <p/\/2) with 
tree-level potential 



Vo{<l>) 



-m 2 <j) 2 + 'j<f)\ 



A 



(31) 



for which the vacuum expectation value of the Higgs is given by v = = 246GeV, 

and A fixes the Higgs mass, m 2 H = 2Xv 2 . Imposing the renormalization conditions that the 
minimum of the potential and the mass of <fi do not change with respect to their tree-level 
values [33], the one- loop zero-temperature correction is given by 



Vi(<f>) 



\- ±9i 

2^ 64vr 2 



m\{(j>) lo 



?» 2 (0) 

m 2 {v) 



2m 2 ((p)m 2 (v) 



(32) 



where gi is the number of d.o.f. of each particle species, m 8 (0) is the ^-dependent mass, 
and the upper and lower signs correspond to bosons and fermions, respectively. The one 
loop finite-temperature corrections to the free energy density, including the resummed 
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Figure 9: The regions of the plane (a c ,«Ar) where detonations and deflagrations may or 
may not exist. The points correspond to the model of section HJ for different values of the 
d.o.f. g and coupling h. Triangles correspond to g = 2, circles to g = 6, and squares to 
9 = 12. 

daisy diagrams are 

^(0,T) = Y^^T J™ dxx 2 log l^exp (-yWm? (0) /T^j 

+ E C ^ 3 (0) " ^ ' (33) 

bosons 

where the upper sign stands for bosons, the lower sign stands for fermions, and M 2 (0) = 
mf (0) + IIj (T), where IT (T) are the thermal masses. The last term receives contributions 
from all the bosonic species except the transverse polarizations of the gauge bosons. 
Hence, the one-loop finite-temperature effective potential is of the form 

F(<f>, T) = V (0) + V x (0) + T) + p A , (34) 

where we have added a constant p\ = — Vo(v) — Vx(v), so that at T = we will not have a 
vacuum energy density of order t> 4 in the true vacuum. Thus, we have p_(T = 0) = 0. If 
all the masses rrii vanish in the symmetric phase, the constant p\ gives the false vacuum 
energy density p+{T = 0). 

For the SM, the relevant contributions come from the Z and W bosons and the top 
quark. The top quark contributes with g t = 12 fermionic d.o.f., with a mass m t (<f)) = 
h t 4>/V2, where h t is the Yukawa coupling (h t /y/2 ~ 0.7). For the weak gauge bosons, it 
can be seen numerically [30] that a good approximation for the one-loop correction is to 
consider a single mass of the form (0) = hb<f>, with hb ~ 0.35, and gb = 6 bosonic 
d.o.f. The rest of the SM d.o.f. have hi <C 1 and only contribute a ^-independent term 
— 7r 2 5 , Hght^ 4 /90, with ^Hght ~ 90. Each extra complex field S gives contributions to the free 
energy of the form (132)330 . with g$ = 2 d.o.f. and a mass of the form m 2 (0) = h 2 (j) 2 + p 2 . 
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The thermal mass is given by II = h 2 T 2 /3 + As/4. For simplicity, in this work we 
shall set \i — \$ = 0. For illustrative purposes, we shall also fix the Higgs mass to the 
value run = 125GeV. It is well known that increasing the Higgs mass weakens the phase 
transition. Therefore, higher masses give lower values of v w . The same happens for nonzero 
values of the boson mass parameter \i since, as /i is increased, the extra bosons decouple 
from the thermal system and the phase transition becomes weaker. In a forthcoming 
paper [35] we will consider several particular (and physically motivated) extensions of the 
S.\I. 

4.2 Critical temperature and latent heat 

We can find the symmetry-breaking minimum (j) m (T) for the free energy ( 1341) by de- 
manding 9jF(0,T)/90 = 0. Then, as we have seen in section dj we define the functions 
JF + (T) = JF(0,T) and JF_(T) = jF(0 m (T),T). From these, we calculate the critical tem- 
perature through the condition JF + (T C ) = JF_(T C ), and the latent heat by 

L = T C (F'_(T C )-F' + {T C )). (35) 

Analytic expressions can be obtained by considering the high-temperature expansion of 
the thermal integrals appearing in (133|) . However, this approximation breaks down for 
rrii((j))/T ^> 1. In our model this corresponds to strong phase transitions, with h<ft/T ^> 1. 
Therefore, we shall calculate numerically T c and L. 

The energy density of the high-temperature phase can be derived from the free energy 
density fl34l) by the relation p + = J 7 + (T)—TJ : ' + {T). The thermal energy density is given by 
p + = p + — pa- In general, we have p + ~ ir 2 g*T A /30, where g* is the number of relativistic 
d.o.f. Thus we can readily calculate the parameter a c = L/p + (T c ). In order to calculate 
the parameter ajy = L/p + (T^), we need to compute the temperature T/v at which bubbles 
nucleate. The latent heat parameter L in oln must be the same as in a c , since for the bag 
equation of state the released energy does not depend on temperature. Therefore, in order 
to apply the results of the previous sections, we must use Eq. ( |35l) for the computation 
of «at, which gives the correct relation a^ja c = (T C /T N ) 4 . If we used instead the energy 
that is released at T = T/v, which is larger, we would be overestimating the velocity, since 
for the bag EOS this would be equivalent to considering a stronger supercooling. 

4.3 The Nucleation temperature 

In a first-order phase transition, the nucleation of bubbles is governed by the three- 
dimensional instanton action 

POO 

S 3 = An / r 2 dr 
Jo 

where V(<fi, T) = JF(0, T) —^(0, T). The bounce solution of this action, which is obtained 
by extremizing 5*3, gives the radial configuration of the nucleated bubble, assumed to be 
spherically symmetric. The action of the bounce coincides with the free energy of a critical 



(36) 



18 



bubble (i.e., a bubble in unstable equilibrium between expansion and contraction). The 
bounce solution obeys the equation 



with boundary conditions 

^(0) =0, lim0(r) = 0. (38) 
The thermal tunneling probability for bubble nucleation per unit volume and time is [9] 

r (T) ~ A (T) e ~ 53(T)/T , (39) 

with A (T) = [S 3 (T) /2vrT] 3/2 . 

The nucleation temperature T/v is defined as that at which the probability of finding 
a bubble in a causal volume is 1, 

dtT (T) V = 1, (40) 



where t c is the time at which the Universe reaches the critical temperature T c and is 
the time at which the first bubbles are nucleated in a causal volume V. In the radiation- 
dominated era we have V ~ (2t) 3 , and the time-temperature relation is given b}{§ 

dT/dt = -HT, (41) 

where H is the expansion rate, H = ^8ttGp + (T)/3. Here, G is Newton's constant. If 
p + w p + ~ 7r 2 g*T 4 /30, then the time-temperature relation is given by the usual expression 
t = £M P /T 2 , where M P is the Planck mass and f = A /45/(167r 3 5'*). 

The nucleation rate T(T) can be calculated by solving numerically Eq. ( 1371) for the 
bubble profile, then integrating Eq. (136]) for the bounce action and, finally, using the 
result in Eq. (|39|) . Analytical approximations to Eqs. (I36H39P have large errors due to 
the exponential dependence of T and the sensitivity of S3 to temperature. We solved Eq. 
( 157)) iteratively by the overshoot-undershoot methocQ. 

4.4 The friction coefficient 

The effect of microphysics on the propagation of the bubble wall can be calculated by 
considering the equation for the Higgs field in the hot plasma. From energy conservation 
considerations, one can derive the equation 



3 This relation may change due to reheating during the development of the phase transition, for t > tpj 

[11301 M- 

4 See Ref. EDI for details. 
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where the sum is over all particle species that couple with 0, m, are the ^-dependent 
masses, and Sfi are the deviations from the equilibrium distribution functions, induced 
by the motion of the wall. The deviations Sfi have been calculated either using kinetic 
theory [TTJ [121 131 EH], or considering infrared excitations of bosonic fields [15], which 
undergo overdamped evolution [T6j . 

Calculations of the friction are usually carried out ignoring hydrodynamics (i.e., disre- 
garding temperature and velocity profiles of the fluid). Assuming stationary motion in the 
z direction, the first term in Eq. (14*2]) becomes (v%, — 1) <p" (z), and the deviations Sfi near 
the wall depend on <p(z) and cj)' (z). To the lowest order, Sfi is proportional to the wall 
velocity v w . As a consequence, if we multiply by <fi' (z) and integrate over — oo < z < +oo, 
we obtain 

^(O,T)-^(0 m ,T)=^ w , (43) 

where 77 is a friction coefficient which depends on the particle content of the plasma. For 
particles with a thermal distribution we have [TT] 

\- Sgjhf f \og Xi \ 2 <p 2 m o- . , 

^ ~ 2. (iyio-iT) \~2^) "T"' (44) 

where Qi is the number of degrees of freedom of species i, hi is the coupling to 0, Tj 
are interaction rates which are typically < 10 _1 T, Xi = 2 for fermions and Xi — f° r 
bosons, while infrared bosons give a contribution [16] 

~ Yl 4r (~^r) l °s(mb((prn) L w ) 1 —, (45) 



where gj, is the number of bosonic degrees of freedom, L w is the width of the bubble wall, 
and mo is the Debye mass, given by m 2 D ~ h\T 2 where h\, is the coupling to <p. 

For bosonic d.o.f., expressions (j44H and fj4"5"]) dominate in different ranges of the cou- 
plings /it, and we shall assume rj = ^ tn + i] iT . Both expressions involve several approxima- 
tions and have 0(1) errors, but should be parametrically correct. 



4.4.1 Parametrization of the friction 

To relate the friction coefficient 77 obtained from microscopical calculations to the one we 
used as a free parameter in the previous sections, compare Eqs. (I43H45I) with Eqs. (|8|9|) 
|36j . We see that the thermal friction coefficient (1441) can be written in the form r] = fjT c o~, 
with 

^ 3gihj ( 'log XiY f (f>m\ 2 (Aa , 

The surface tension is roughly given by a ~ <p 2 rn jL w . Furthermore, in general we have 
4> m ~ T. Therefore, the last factor in Eq. f|4"5]) becomes T 3 /L w ~ Ta. Thus, the infrared 
boson contribution can also be written in the form (Q, with 




gauge 



Thus, for a given model, we can estimate the value of 77 in Eq. ([9]) using Eqs. (I46|47l) . 
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4.5 Numerical results 



We have numerically calculated the values of a c , a at and T]/L for this model, as functions 
of the coupling h and the number of d.o.f. g of the scalar singlets. From these, we have 
calculated the wall velocity v w at T = T N , using the results of section [31 

In Fig. [10] we have plotted the value of the deflagration wall velocity as a function of 
h for g = 2. Since the phase transition strengthens with h, one would expect the wall 
velocity to be a monotonically increasing function of h (perhaps eventually bounded by 
c g ). Indeed, roughly we have v w (Tat) ~ AjF(T/v) /rj, so v w is dominated by the amount 
of supercooling and by the friction. The friction coefficients ??th and r] iT oscillate for h ~ 1 
due to the log terms in Eqs. (pHj) and fT45|) . This causes a minimum and a maximum in 
the wall velocity. For large h, the friction coefficient r^ t h dominates since it goes like h 4 , 
and the velocity finally decreases. To illustrate this effect, we have plotted in Fig. fTOlthe 
value of the friction coefficient f] = fj t h + fjh- 
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Figure 10: The wall velocity and friction coefficient as functions of the coupling h for 
g = 2 d.o.f. 

Figure [TTJ shows the deflagration wall velocity as a function of the coupling h for the 
cases g = 2, g = 6, and g = 12. The crosses in the g = 12 curve indicate the points 
where <p (T/v) /TV ~ 1 and (T/v) /TJv ~ 2. (These points were taken as a reference for 
the variation of the parameters in the previous section.) For each g, the calculation was 
done up to a value h = /i max for which the time required to get out of the supercooling 
stage becomes too long for the numerical computation. This happens because the phase 
transition becomes so strong that the barrier between minima persists at low tempera- 
tures, and the nucleation temperature quickly falls to zero. Indeed, near this endpoint, 
the free energy has a barrier between minima already at T = 0. For h beyond h = h m3/X , 
the temperature T c also goes to 0. The dependence on the coupling will not be so strong 
if fermions are added to the model (see e.g. Fig. 6 of Ref. [30]). 

Notice that the deflagration wall velocity can reach the speed of sound c s only for 
large enough g and the largest values of h in each curve, corresponding to really strong 
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Figure 11: The wall velocity function of h. 

phase transitions. Immediately after crossing the line of v w = c s , the deflagration solution 
may disappear (cases g = 6 and g = 12 in the figure), and it may reappear (case g — 6), 
indicating a borderline case in parameter space. For the case m# = 125GeV, we have 
considered values of in the range — 20, varying h up to /i max for each g, and we have 
not found detonation solutions. This occurs because, although the supercooling increases 
with h, the friction increases as well. For lower Higgs masses, namely, m# = lOOGeV, 
detonation solutions appear only for h « h max . 

Figure M shows the values of a c and corresponding to the curves of Fig. [TTJ Lower 
values of h give weaker phase transitions and, hence, lower values of a c and a^. As h 
becomes large, the supercooling diverges before a c reaches the limiting value 1/3. Most 
points, though, are close to the region where detonations do not exist. For these cases, 
detonations would exist only for very low values of the friction. On the other hand, 
for large h the points move away from the lowest curve in Fig. [91 For these points, 
the supercooling is quite larger than the minimum that is needed for the existence of 
detonations. Therefore, one expects that this kind of solution would appear if the friction 
were a little smaller. As we already mentioned, the error in the calculation of the friction 
is a (9(1) factor. Thus, we also consider a friction which is a factor of 2 larger and a factor 
of 2 smaller than the value given by Eqs. (I44II45I) . The result is shown in Fig. [T2J In the 
case rj/2 we obtain detonations for large values of h, if g is large enough. 

5 Analytic approximations for the wall velocity 

Analytic approximations for the wall velocity are very useful for cosmological applications. 
Two approximations are commonly used. In the case of deflagrations, if one ignores 
hydrodynamics, the wall velocity can be obtained by equating the pressure difference 
Ap between phases (which pushes the wall towards the high-temperature phase) to a 
friction force per unit area, which is typically proportional to the velocity, / = r\v w . 
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Figure 12: Same as for Fig. [TTJ but considering a larger and a smaller friction parameter 
(from left to right, the curves correspond to g = 12, 6, and 2). 



Thus, one obtains the well known result v w = Ap(T)/rj, which has been widely used in 
electroweak baryogenesis [H] [TJl EH HU EES, E32S EH] • This velocity depends on the amount 
of supercooling, since Ap = at the critical temperature, and gives a good approximation 
for nonrelativistic velocities. However, ignoring hydrodynamics overestimates the value 
of v w , since the released latent heat accumulates in front of the wall, causing a slow-down 
of the wall velocity. Nevertheless, if the friction is large enough, this effect is negligible. 

In the case of detonations, the Chapman- Jouguet condition, t>_ = c s , leads to the 
simple formula 

vg3 + Vg+2g3 
1 + a 

which is obtained by setting t>_ = in Eq. This formula has been widely used 

for the calculation of gravitational radiation [SI El EH- However, as we have already men- 
tioned, it is rather strange that Eq. (1481) does not depend on microphysics. Furthermore, 
it depends only on the ratio L/p + (Tjy). This means that increasing the latent heat will 
cause the same effect as decreasing T N , which is to increase the wall velocity It is true 
that the pressure difference goes like the product L(T C — T N ), but, on the other hand, 
increasing L will enlarge the effect of hydrodynamics, which is to slow down the wall. 
As we shall see, the velocity decreases if L is increased with the amount of supercooling 
held fixed 0. As can be seen in Fig. HJ according to our results, the Jouguet velocity is 
a bad approximation to the detonation velocity (it approximates well the non-physical 
solution). 

In this section we shall derive analytical approximations both for the detonation and 
the deflagration, by taking the ultrarelativistic and nonrelativistic limits, respectively, of 
the equations obtained in section [2J Before going on, though, one further comment is 



5 In a particular model, increasing L will give in general more supercooling. However, such a relation 
cannot be derived from hydrodynamics alone, but from calculating the nucleation rate. 
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worth. As we have seen, the parameter e in a = e/p + is not given by the latent heat L, 
but by e = L/4. This fact was already mentioned in earlier works on gravitational waves 
[5]. However, the vacuum energy density e is sometimes confused with the latent heat 
L. When analytical results are used for a particular model, this amounts to considering 
a value of a four times larger than the actual value, which leads to an overestimation of 
the detonation velocity. 



5.1 The ultrarelativistic limit 

For v w ~ 1, we can write v w — 1 — 8, with 5 <C 1. In this limit we must consider 
detonations, for which v w = \v + \ and a + = = ot. Therefore, we have \v + \ = 1 — 5 and, 
from Eq. (jHJ), = 1 — 5 (1 + 3a) + (5 2 ). To lowest order in 5, we have |f+| 7+ = 1/V25 
and |w_|7_ = 1/^25 (1 + 3a). Thus, the last term in Eq (JHD is ~ 1/V6. For the first 
term we have, from Eq. f[2"Tj) . (p + — p_) / p+ = —2a + O (5). The second term is obtained 
from Eqs. (1221) and (j23|) . From Eq. (|23"1) we have p~/p+ = 1 + 3a + O (5), from which we 
obtain T_/T + = (a_/a + )" 1/4 (1 + 3a) 1/4 + O (5), and = (a/a + ) 1/4 (1 + 3a) 3/4 + 

O (5). Inserting these results in Eq. (|20|) we obtain, to lowest ordei|j, 

v w = 1-5, 

26 = (3/4)(yTTto + l)^ 

r 1 /4(l + 3 «) 5 / 4 _ r -i/4(i + 3a )3/4' 

with a = e/ (a+T^-), fj = r]/ (a + T^), and r = a„/a + = 1 — 3a c , with a c = e/ (a + T c 4 ). We 
remind also that the correct relation between £ and L is e = L/4. 

Since the detonation velocity is usually large (v w > 0.8), the simple formula (149]) turns 
out to be an excellent approximation in general. We have plotted this approximation in 
Fig. [13] together with the numerical result, for the same detonation curves of the right 
panel of Fig. [6] For comparison, the Jouguet velocity lies below 0.7 (and is constant for 
a fixed). Since the parameter 6 is bounded by 1 — c s ~ 0.4, we expect Eq. (T49]) to give a 
good approximation in all the possible range of the detonation velocity. 

It is interesting to consider some limiting cases of Eq. (T4*9]) . For rj — > with all the 
other quantities held fixed, the wall velocity approaches the speed of light. However, for 
any finite value of rj, the denominator diverges at / 'a c = a + / 'a_ = 1/ (1 — 3a c ). Hence, 
a 7v has a lower bound larger than a c (except in the limit a c — > 0). This divergence reflects 
the fact that a certain amount of supercooling is always needed for detonations to exist. 
As a consequence, we cannot take the limit of small supercooling, a^ja c ~ 1. 

If we increase L with fixed qjjv and r], it can be easily seen from Eq. ( 149]) that v w 
decreases. This is because increasing L strengthens the effect of hydrodynamics, which is 
to slow-down the wall. Taking the limit of large L is not very useful. In the first place, L 



6 Keeping the next order in the equation is straightforward, and leads to 

yJfD* - f] 2 (Aq + 3gV2 + l _ 2/q)(q^ + 1) - 

(7j/2)(4g + 3gV2 + i_2/ 9 ) ' 

where D is the denominator in Eq. (|49[) and q = 1 + 3a. 
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Figure 13: The analytical approximation (|49|) (dashed line) together with the numerical 
result (solid line), for the detonation solutions of Fig. El right. 

is bounded by (4/3)p + (T c ) (we remark that this is a general thermodynamical constraint 
|30j). Hence, this limit corresponds to taking a c — > 1/3 (or a_ — > 0). However, this is 
the limit of an extremely strong phase transition which, in a realistic model will have also 
Tjv — > (or «at — > oo). Conversely, the limit of large supercooling, a^ja c — > oo, should 
be taken together with a c — > 1/3. It is not straightforward how to take this limit in Eq. 
( H9|) . since the bag model does not provide a way to calculate the nucleation temperature 
T N . 

On the other hand, we can consider the limit L — > 0. In this case both a c and — > 0, 
(unless T/v <C T c , which will not be the case in a realistic model if L is small). Hence, 
Eq. (@SD yields « (p/p+) / ( ajv - a c ) = (Aij/L) / (1 - T 4 /T 6 4 ). From Eqs. (HKUD, 
we see that L (1 — T^/T c 4 ) /4 = p_ (T N ) — p + (T N ). Since in the ultrarelativistic limit we 
have l/\/25 = v w j w , we obtain v w ^ w = Ap(T N )/rj. Thus, as expected, in the limit L — > 
we can neglect the effects of hydrodynamics. The wall velocity becomes, according to Eq. 
P}, v w = l-ifa/ Ap(T N )f. 

5.2 The nonrelativistic limit 

In the nonrelativistic limit ^ < 1, we must consider deflagrations. Thus, we have 
\v~\ = v w , and Eq. (j3J) becomes i> + = (l — 3a + )v_ + 0(v^). Therefore, according to 
Eq. (J2U), we have (p+ - p_) /p + = O (v 2 J and, by Eq. (J23D, p_/p + = 1 - 3a + + O ). 
Thus, T_/T + = (a_/a + )~ 1/4 (1 - 3a + ) 1/4 + £> (v*) an d s_/s + = (a_/a + ) 1/4 (1 - 3a+) 3/4 + 
C (v^). To lowest order in v w , Eq. fTSOj) becomes 

x X 1/4 xx -1/4 

3a + (2 - 3a+) = f (1 - 3a + ) 3/4 - f (1 - 3a+) 1/4 + 3a+. (50) 

In the case of deflagrations, we must also find function of a N . In this approx- 

imation, v + — V- is given by 3a + v w . Hence, according to Eq. ([30]) . a a? — a + is of order 
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v w , and so the left hand side becomes \f?> (a/v — a+) /(4ajv). Thus, to lowest order in t>„,, 
Eq. ( 1301) gives a + = «at — 4v^3a 2 v tv Inserting this in Eq. ([50]) and discarding higher 
order terms, we obtain 

v w = 

(1 - 3a)^ ^3a - r~i (1 - 3a)^ + A (1 - 3a)*) 
fry (2 - 3a) (1 - 3a)* + Sv^a 2 

where a = a^ = (L/4) / (a + T^), fj — r]/ (a+Tjy), and r = a_/a + = 1 — 3a c , with 
a c = (L/4)/(a + T 6 4 ). 

The appearance of factors (1 — 3a7v) 1//4 indicates the fact that the deflagration solution 
cannot exist for a at arbitrarily large. For r] — > 0, only the first term in the denominator 
vanishes; the velocity will not necessarily be large, nor the approximation break-down. 
This reflects the fact that hydrodynamics slows down the motion of the bubble wall. For 
L — > 0, i.e., a c ,a7v — > 0, we obtain, similarly to the detonation case, v w = Ap(T^)/r], 
which corresponds to neglecting hydrodynamics. It is interesting to consider also, without 
neglecting hydrodynamics, the case in which Ap is small, i.e, the limit of small supercool- 
ing. 

5.2.1 Small supercooling 

In the nonrelativistic limit, the pressure difference p + — p- is O (v^). This can be seen 
already from the nonrelativistic version of Eqs. (pQ), 

w + v + = W-V-, w + v^_ + p+ = w^v 2 _ +p-. (52) 

On the other hand, according to Eq. fflBl . the temperature difference T + — T_ is O (r]V w ). 
We shall now consider the case in which rj is not too small, so that the temperature 
difference must be T + — T_ = O (v w ). This means that, to order v w , the pressure difference 
due to supercooling, p_ — p + ~ T c — T, is canceled by the pressure difference due to the 
temperature gradient, p- — p + ~ T_ — T + . Hence, the differences T c — T± must be also 
O (v w ). Thus, in this case we can expand the thermodynamic quantities to first order in 
T—T c . In particular, for the pressure functions we have p± (T) = p± (T c )+s± (T c ) (T — T c ), 
which can be used instead of the equations of state (FIB]) . Thus, we can write 

p + (T+) - p- = (s + - s_) (T + - T c ) + S _ (T + - T_) , (53) 

where s + and s_ are evaluated at T = T c . Therefore we have 

T + -T_ = ^(T c -T + ) + 0{vl). (54) 

To first order in v w , we neglect the first term in Eq. (TlBl and, inserting Eq. (|54l) in the 
second term, we obtain 

T c -T + s_ 

L — ^ = —r]v w , (55) 

T c s + 



4 +4(1 -3a)* -3r3 (1 - 3a 
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- w • 



where we have used L = T c (s+ — s_) and v + = (w-/w + )v~, with = v v 

If we had T + = Tn, this would only give a correction of a factor to the usual 

equation v w = Ap (T) / 7]. However, an important contribution comes from the fluid that 
is accumulated in front of the moving wall. For T + and ph T c , the relation (130]) between 
the nucleation temperature T N and that of the reheated fluid T + is given by 

T+ — T N _ v w L 
T c \/ r 3w + 

where the enthalpy w + is evaluated at T = T c . We finally obtain 

(w+/wJ) (T c - T N ) /T c 



T]/L + L/{y/ZwJ) 



(57) 



which agrees with the result of Ref. [10] . Comparing with Eq. flSj) and taking into account 
that AjF (T) pa L (T c — T) /T c , we may write 

L (T c - T N ) /T e , \ 

v w = —-^ 58 

^cff 

where the effective friction coefficient 

Vez = -(v + -^-) (59) 



includes the effects of hydrodynamics. If friction dominates, i.e., if rj/L ^> L/(v / 3u'_), 
we obtain v w = (w-./w+)L (T c — Tjv) /^. For L p + , we have w_ ~ w + , and we recover 
Eq. (JH1). In the opposite limit, i.e., for rj/L <^ L/ (\/3w_), we still obtain a finite velocity 

= (4p + /L) (T c — T N ). Again, this means that the latent heat that is accumulated in 
front of the bubble wall slows down the motion of the wall. 

In Fig. [14] we plot the two approximations for deflagrations, together with the nu- 
merical result and the usual approximation which neglects hydrodynamics. As expected, 
the numerical calculation gives smaller values of v w than the analytical approximations, 
since the exact velocity has an upper bound, namely, the speed of light. We see that in 
general the nonrelativistic approximation is quite good up to v w ~ 0.4, and breaks down 
for v w > c s . Notice that a useful, rough approximation consists of assuming that v w is 
given by, say, Eq. f|5T|) up to v w = c s , and by v w = c s when Eq. flBTI) becomes supersonic. 



6 Conclusions 

We have studied the steady state velocity of bubble walls in first-order cosmological phase 
transitions. We have numerically investigated the whole velocity range < v w < 1, and 
we have found analytic approximations for the limits t) u <l and v w ~ 1. 

Thus, we have considered both detonations and deflagrations. We have derived a 
relatively simple set of equations for the wall velocity and the fluid variables on both 
sides of the wall. To obtain these equations we have used three ingredients, namely, the 
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Figure 14: The various analytical approximations for deflagrations. Solid line: the numer- 
ical result. Dashed line: the nonrelativistic approximation of Eq. (I5T1) . Dashed-dotted 
line: the approximation of Eq. (1571) for small supercooling. Dotted line: the usual ap- 
proximation v w = Ap/rj. 

continuity conditions for the fluid variables at the phase transition front (and at the shock 
front in the case of deflagrations), the introduction of a damping term in the equation 
of motion for the order parameter (Higgs field), and the bag equation of state. The 
parametrization of the friction provides a model-independent equation which involves a 
friction parameter rj. We obtained this equation by generalizing a procedure used in Ref. 
[TO] for the case of nonrelativistic velocities and small supercooling. We have also shown 
how to calculate the parameter rj in a specific model. 

We have solved numerically these equations for the wall velocity as a function of the 
parameters a c , and rj/L which depend essentially on the latent heat L, the nucleation 
temperature T/v, and the friction coefficient rj. We have studied the regions in parameter 
space where detonations and deflagrations are possible, and we have also considered a 
specific model (the Standard Model with extra singlet scalar fields) to investigate how a 
physical case moves about through these regions. 

We have found that supersonic velocities (either detonations or deflagrations) are in 
general obtained only for very strong phase transitions, with (T/v) /T/v > 2, occurring in 
the model for relatively large values of the number of scalar singlets g and of their coupling 
h to the Higgs. Besides, the existence of detonations requires relatively low values of the 
friction, and it may happen that for such strong phase transitions neither detonations 
nor deflagrations exist. In such a case, the steady state cannot be reached, and the wall 
will accelerate until bubbles collide. The case in which the wall velocity gets close to the 
speed of light before bubble collision may have important implications for gravitational 
wave generation. Interestingly, this ultrarelativistic situation has been considered very 
recently [37] for the SM extension with singlet scalar fields, finding that such "runaway" 
solutions exist for very strong phase transitions. 

In a more general extension of the SM, the singlets may have a ^-independent mass 
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term, which will weaken the phase transition. In that case, detonations may not appear 
at all. The inclusion of fermions in the model will also weaken the phase transition, but 
the friction coefficient will not have the infrared contribution and will be smaller. In a 
forthcoming paper [35] we shall investigate several extensions of the SM, including the 
MSSM and an extension with strongly coupled fermions [38J. 

Finally, we have discussed the validity of commonly used analytical approximations 
for the wall velocity, and we have obtained nonrelativistic and ultrarelativistic approxi- 
mations. In particular, we have shown that the Jouguet velocity gives a really bad ap- 
proximation to detonations; the actual detonation velocity is quite larger. We have found 
an alternative approximation (j4"9"|) which, besides being more realistic, is quite simple and 
gives an excellent fit to the numerical result. 
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